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We present a comparative full-potential study of generalized Kohn-Sham schemes (gKS) with explicit focus 
on their suitability as starting point for the solution of the quasiparticle equation. We compare GqWo quasipar- 
ticle band structures calculated upon LDA, sX, HSE03, PBEO, and HF functionals for exchange and correlation 
(XC) for Si, InN and ZnO. Furthermore, the HSE03 functional is studied and compared to the GGA for 15 
non-metallic materials for its use as a starting point in the calculation of quasiparticle excitation energies. For 
this case, also the effects of selfconsistency in the GW self-energy are analysed. It is shown that the use of a 
gKS scheme as a starting point for a perturbative QP correction can improve upon the deficiencies found for 
LDA or GGA staring points for compounds with shallow d bands. For these solids, the order of the valence and 
conduction bands is often inverted using local or semi-local approximations for XC, which makes perturbative 
Go Wo calculations unreliable. The use of a gKS starting point allows for the calculation of fairly accurate band 
gaps even in these difficult cases, and generally single-shot GqWo calculations following calculations using the 
HSE03 functional are very close to experiment. 

PACS numbers: 71.15.Mb, 71.15.Qe, 71.20.Nr 



I. INTRODUCTION 

Density functional theory (DFT) has become the most suc- 
cessful method for condensed matter calculations. This suc- 
cess is largely rooted in the simplicity of the exchange and 
correlation (XC) energy in the local density (LDA) or gener- 
alized gradient (GGA) approximation. However, the underly- 
ing Kohn-Sham (KS) formalism fails in the prediction of elec- 
tronic excitation energies of semiconductors and insulators.' 
A significant step forward to correct excitation energies was 
achieved when the first ab initio calculations of quasiparticle 
(QP) states were performed. ^'^ Their description is based on 
the quasiparticle equation with the XC self-energy E for one 
excited electron or hole i''^'^'^ In general, its solution is based 
on Hedin's GW approximation (GWA) for the self-energy^ 
and a perturbative treatment of the difference to the XC po- 
tential used in the KS equation. The central quantity is the 
dynamically screened Coulomb potential W, which character- 
izes the reaction of the electronic system after excitation. For 
many non-metals, such as the semiconductor silicon (Si), the 
method works well with an accuracy of 0.1-0.3 eV for their 
QP gaps,' if the Green's function G is described by one pole 
at the KS (i.e.. Go), or better, at the QP (i.e., G) energy. How- 
ever, for systems with a wrong energetic ordering of the KS 
bands first-order perturbation theory is not applicablci^ Ex- 
amples are semiconductors with a negative fundamental gap 
in DFT-LDA or -GGA, e.g. InN^ (and references therein), or 
with shallow d bands, e.g. ZnOii 

Besides the KS approach itself, the origin of the band 
gap problem is related to the semi-local approximation 
(LDA/GGA) for XC, which introduces an unphysical self- 
interaction and lacks a derivative discontinuity.^''^ These defi- 
ciencies can be partially overcome using self-interaction-free 
exact exchange (EXX) potentials,'" which are special realiza- 
tions of an optimized effective potential (OEP) method.^ A, 



by conception, different way to address the band gap problem 
is the use of a generalized Kohn-Sham (gKS) scheme, which 
means starting from a scheme with a spatially non-local XC 
potential." In this framework, the screened-exchange (sX) 
approximation uses a statically screened Coulomb kernel in- 
stead of the bare kernel in the Hartree-Fock (HF) exchange— 
and, with it, resembles the screened-exchange (SEX) con- 
tribution to the XC self-energy in the GWA.^"* Other hy- 
brid functionals such as those following the suggestions of 
Adamo and Barone (PBEO)'^ or Heyd, Scuseria, and Ernzer- 
hof (HSE03)'^ combine parts of bare or screened exchange 
with an explicit density functional. The gKS eigenvalues are 
usually in much better agreement with the experiment than 
the LDA/GGA ones. Therefore, the gKS solutions are sup- 
posed to be superior starting points for a QP correction, since 
first-order perturbation theory should be justified. Hence the 
replacement of G by Go calculated from solutions of a gKS 
scheme may be interpreted as a first step towards a self- 
consistent determination of the self-energy operato r '^''^ 



Here we report a systematic study of QP energies calculated 
from GW corrections to the results of gKS schemes. First, in 
section IIIIAI we evaluate the performance of GqWo correc- 
tions to different gKS starting points for Si, InN, and ZnO 
using the functionals sX, HSE03, PBEO, and HF. The results 
are compared to those of the standard KS approach based on 
an LDA functional. In section llllBI QP gaps are calculated for 
a benchmark set of fifteen non-metalls utilizing the HSE03 
starting point. The effects of selfconsistency in G and W are 
discussed in comparison to the results based on a GGA start- 
ing point used in Ref. il6i. 



II. METHOD 

All calculations are performed at the experimental lattice 
constants. We use the projector augmented-wave (PAW) 
method as implemented in the Vienna Ab initio Simulation 
Package. ^^ For the details of the GW implementation we re- 
fer to Refs. [181 and 116|. The GW calculations are carried out 
using a total number of 150 bands for all materials. For the 
Brillouin-zone integrations, 8x8x8 k-point meshes includ- 
ing the r point were used, except in the case of ZnO (LDA) 
and InN (sX, HSE03) where the k-point convergence of W 
was found to be critical for meshes containing P. In these 
cases, 8x8x8 Monkhorst-Pack k-point grids avoiding F 
were used for the evaluation of W. 

One problem of the presence of shallow d levels is the 
strong core-valence XC interaction.''*''^ It can be estimated 
within the LDA or HF approximation, where the latter one 
is expected to be more reliable since the GW self-energy ap- 
proaches the bare Fock exchange operator in the short wave- 
length regime (i.e., at large electron binding energies). There- 
fore, we apply the HF approximation to the core-valence XC 
self-energy for all GW calculations reported here (see Ref. il6i 
for a validation of this approach). 

Here we do not update the QP wave functions correspond- 
ing to the neglect of non-diagonal matrix elements of the self- 
energy represented in terms of the gKS wave functions i//| . 
The QP excitation energy £^+' of a state A in the (NH-l)-th 
iteration is related the N-th iteration through the Unearized 
equation: 
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where T is the kinetic energy operator, V„-e the nuclei- 
potential, Vh the Hartree potential, and Z^ the renormaliza- 
tion factor given by 
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The first iteration, usually denoted by GqWo,'''^ is based on 
the gKS eigenvalues ej = ef as initial input to the GW cal- 
culation. Within this approximation the perturbation operator 
in ([TJ becomes E(ef ) — VJc corresponding to the differ- 
ence between the GW self-energy and the non-local XC po- 
tential used in the gKS equation, and the GqWo-QP shift for a 
certain gKS state is given by: 



potential. The superscript DFT indicates that the respec- 
tive quantity is evaluated in some (quasi)local approximation, 
while E^'' corresponds to one of the (screened) Coulomb ker- 
nels given in Table H] The weight a of the short range part and 
the inverse screening length ji are also listed in this table. For 
simplicity, in the case of sX, the inverse screening length jj. 
corresponding to the Thomas-Fermi wave vector was chosen 
materials independent kjF — 1-55 A^'. Gradient corrections 
were used for the HSE03 and PBEO functionals. 



III. RESULTS 
A. Go Wo quaslpartlcle band structure of Si, ZnO, and InN 

1. Generalized Kohn-Sltam bands 

Results of the gKS calculations are summarized in Table Ull 
together with experimental results. We notice an increase of 
the computed gaps when going from the LDA to a truly non- 
local XC functional. This can be attributed to the reduction 
of the spurious self-interaction found in LDA, the inclusion 
of a potential discontinuity between filled and empty states in 
the XC functional j^'^'" and the - in comparison to LDA - en- 
hanced core-valence exchange. Furthermore, for direct gaps 
we observe the tendency to increase from LDA over HSE03, 
sX, and PBEO to the HF values. 

For ZnO and InN, the tZ-band binding energies increase 
with respect to the LDA, approaching the experimental val- 
ues. This is an important fact to note, since the difficulties 
of the LDAh-GoWq approach for these compounds partly re- 
sult from too shallow d electrons in LDA, which is itself a 
result of the spuriously contained self-interaction in local and 
semi-local functionals. The too shallow d electrons hybridize 
too strongly with the p bands at the valence-band maximum 
(VBM), pushing them upwards, in turn decreasing the gap (pd 
repulsion)'' beyond the common LDA gap underestimation 
found for example in Si. Contrary to the LDA, the gKS start- 
ing points yield stronger bound d electrons (cf. Table HJi with 
binding energies closer to the experimental values. Hence the 
gKS functionals can be expected to give a more reasonable 
estimate for the influence of the pd repulsion on the funda- 
mental gap and to provide better starting wave functions for 
the GW calculations. For instance, in ZnO the Zn d character 
of the wave functions at the VBM is decreased from 0.3 for 
LDAto0.25forHSE03. 

In comparison to experiment the gKS functionals, with ex- 
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In the actual implementation of the gKS schemes, we split 
the gKS XC energy into the form: 
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i.e., a short-range non-local exchange term is added and 
treated exactly resulting in a non-local (screened) exchange 



TABLE I: Parameters of the gKS exchange functionals used. 



Functional 



.ST Coulomb kernel 



M[A-'] 



LDA 


0.00 


l/|x| 


- 


sX 


1.00 


exp(-M|x|)/|x| 


L55 


HSE03 


0.25 


erfc(/i|x|)/|x| 


0.3 


PBEO 


0.25 


l/|x| 


- 


HF 


1.00 


l/|x| 


- 



TABLE II: Direct and indirect generalized KS band gaps E",, .. and 

average d-band binding energies E^ calculated for cubic Si, InN, 
and ZnO. Experimental values are from data collections in Refs.ja 
116.20 , 21 , 22,23 . In the case of InN and ZnO they refer to the wurtzite 
polytype. All values are given in eV. 





Energy 


LDA 


sX 


HSE03 


PBEO 


HE 


Exp. 


Si 


pgKS 


0.51 
2.53 


0.98 

3.23 


1.04 
3.15 


1.85 
3.98 


6.57 
9.08 


1.17 
3.40 


InN 


pgKS 
pgKS 


-0.38 
13.1 


0.39 
17.2 


0.37 
14.6 


1.14 
14.7 


7.15 
18.6 


0.61 
16.0-16.9 


ZnO 


Er 
Efr 


0.6 
4.6 


2.97 
8.2 


2.11 

5.7 


3.03 
5.8 


11.07 
9.3 


3.44 
7.5-8.8 



FIG. 1: (Color online) First-order quasiparticle shifts A^^^ versus 
gKS eigenvalues for InN (a) and ZnO (b). The valence-band max- 
imum (VBM) and the conduction-band minimum (CBM) are taken 
as energy zeros for occupied and empty states, respectively. Results 
for five different starting gKS band structures are shown. 
a) 3 
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correlation self-energy T.c scaled with the dynamical renor- 
malization factor and consequently undergo the sign change 
of Ec ^t the Fermi level. For the LDA exactly the opposite 
sign change is observed, i.e. valence bands aquire a negative 
shift and conduction bands are shifted upwards. It is remark- 
able that the two hybrid functionals HSE03 and PBEO, which 
essentially mix 25 % HF and 75 % DFT exchange, yield very 
small QP shifts. This highlights that the 1/4 recipe is indeed a 
remarkable good and robust choice not only for total energies 
but also for one-electron energies in semiconductors. 

The performance of GW upon sX is somewhat disappoint- 
ing. The li-band positions were very well described using the 
sX approach (see Table Hill, but applying GW corrections, the 
localized d states are shifted upwards by 1.5-2 eV, signifi- 
cantly deteriorating agreement with experiment. We will re- 
turn to this issue in the next section and in the conclusions. 
Since the SEX term in the GW approximation is very simi- 
lar to the sX term in gKS schemes, the differences must be 
mainly related to a different behavior of the Coulomb hole 
(COH) term in the GW approximation and the local-density 
part of the sX functional. 

Turning now to a detailed discussion of the InN results [Fig. 
[Tta)], we note that the QP shifts for the upper valence bands 
and the conduction bands show a rather weak dispersion with 
the eigenstate for all gKS schemes. In this energy range, they 
may be approximated by scissors operators for occupied and 
empty states, whose absolute values depend on the used gKS 
functional. Only for the lowest In 4d- and N 2s-derived bands, 
at gKS energies below —10 eV, noticeable deviations occur in 
the form of a down or upward shift for the LDA and sX start- 
ing points, respective. In the considered energy range the QP 
shifts calculated for the HSE03 or PBEO starting point are re- 
markably small with respect to the dispersion and amplitude, 
and they remain well below 1 eV. For ZnO [Fig.lTJb)] similar 
observations can be made, however, as discussed before, the 
smaller binding energy of the Zn 3d states causes a stronger 
hybridization of /:> and d states than in InN. This might explain 
the upward bending of the sX shifts close to the VBM and the 
stronger bending of the LDA conduction band shifts. 



■gKS VBM ' 



^gKS-'=CBM (eV) 



ception of HF, generally perform better than standard DFT- 
LDA for the gaps and c/-electron binding energies. For InN 
there is another important fact to note. All gKS functionals 
are found to yield the correct ordering of the Fi^ and Fis,, 
states at the zone center, in contrast to LDA findings which 
give a negative sp gapi^ This is essential for the QP descrip- 
tion, since a correct energetic ordering of the single-particle 
states is an inevitable prerequisite for a perturbative treatment 
of the GW correctionsJii^ 



2. Quasiparticle shifts 

Figure [T] shows the calculated QP shifts [Eq. ^} for InN 
and ZnO plotted versus the gKS eigenvalues. In the extreme 
limit of the HF starting point the QP shifts are given by the 



3. Non-selfconsistent QP bands 

Table Uni shows the results of the GqWq calculations for the 
fundamental gaps and li-electron binding energies. In contrast 
to the LDA or gKS energies, these QP energies should have 
a physical meaning as measurable quantities and thus can be 
compared directly to the experimental gap values. However, 
it is necessary to note that the experimental values for InN 
and ZnO correspond to the wurtzite instead of the zincblende 
polytype used in the calculations. The zincblende gaps are ex- 
pected to be about 0.2 eV smaller than the wurtzite ones.^' — 
For silicon, the QP gaps calculated upon the LDA, sX, or 
HSE03 starting point bracket the experimental values closely, 
with a maximum deviation of 0.2 e V. The LDAh-GqWo values 
slightly underestimate the gaps. In fact, this corresponds to the 
general trend of gap underestimation in the LDA/GGAh-GqWo 
approach using full-potential methods recently established by 



TABLE III: Direct (d) and indirect (i) GqWq QP band gaps, average 
t/-band binding energies and static electronic macroscopic dielectric 
constants £00 calculated upon the respective gKS band structures for 
Si, ZnO, and InN. All energy values are given in eV. Experimental 
results are given for comparison. 





Energy 


LDA 


sX 


HSE03 


PBEO 


HP 


Exp. 


Si 


si 


1.08 


1.31 


1.32 


1.65 


2.93 


1.17 




pQP 


3.18 


3.49 


3.48 


3.72 


5.21 


3.40 




£00 


13.9 


10.8 


9.8 


7.8 


3.4 


11.90 


InN 




0.00 


0.55 


0.47 


0.78 


2.56 


0.61 




pQP 


15.1 


15.6 


15.2 


15.3 


16.6 


16.0-16.9 




eco 


12.2 


6.6 


6.8 


4.9 


2.4 


7.96 


ZnO 


pQP 
^g 


2.14 


3.36 


2.87 


3.24 


5.71 


3.44 




^QP 

d 


5.6 


6.2 


6.1 


6.2 


7.0 


7.5-8.8 




£00 


5.3 


3.0 


3.4 


3.0 


1.8 


3.74 



different groups. '^'^^ For the gKS starting points, the calcu- 
lated gaps are larger than the experimental ones. In more de- 
tail, starting from sX or HSE03, we find virtually identical 
results very close to the experiment, while the PBEOh-GqWo 
approach overestimates the gaps significantly. For the HF 
starting point, the largest deviations are found, which could 
be expected already from the large differences between the 
QP eigenvalues e^ and the HF ones e^^ and the consequent 
break-down of the perturbative GqWq treatment in equ. ([TJ. 
For the PBEO starting point, however, the deviations demand 
for a more subtle interpretation. To understand this devia- 
tion, the static RPA dielectric constants calculated from the 
gKS eigenvalues and wave functions are included in Table Hill 
Indeed, a comparison with the experimental values suggests 
significant underscreening for the PBEO starting point, which 
certainly contributes to the overestimation of the gaps. The 
values calculated upon the sX and HSE03 functional are in 
much better agreement with the experiment, which is brack- 
eted closely by them (sX/HSE03) and the LDA value. 

For InN and ZnO, the trends just discussed for Si still hold, 
but the actual benefit of a gKS starting point becomes more 
apparent. In contrast to the LDA, the gKS starting points re- 
produce the coiTect ordering of gap states for InN and yield 
more meaningful dielectric constants in the case of HSE03 
and, with some restrictions, in the case of sX. Taking into 
account that different polytypes are compared the agreement 
between the experimental gap values and the sX+GqWo or 
HSE03h-GoWo results is fair 

Similar to the gaps, the QP t/-band binding energies calcu- 
lated for InN and ZnO show a much weaker variation with 
the starting XC functional than the original gKS one-electron 
t/-band energies. The agreement of the calculated energies 
with experimental values (Table Hllll is good for InN. For ZnO 
with more shallow d levels the picture is less clear, although, 
here the comparison with experiment is also hampered by the 
remarkable scatter in the measured values. In general, we 
observe that the d states become more strongly bound for 
the gKSn-GoWo approach than for LDAh-GqWo, therefore cer- 
tainly moving in the right direction. Hence, the better treat- 



TABLE IV: Results for the fundamental gaps of the HSE03 and 
quasiparticle (GoWq. GWq and GW) calculations, and static elec- 
tronic macroscopic dielectric constants as used in Wq (RPA). The cal- 
culated values for the spin-orbit coupling (SO) induced gap-closing 
given in the last column have been included in the gaps. Also re- 
ported is the mean absolute relative error (MARE) and the mean rel- 
ative error (MRE) for the gaps. Experimental data for the gaps and 
dielectric constants are given for comparison (for references see Ref. 
16), underlined values indicate zero temperature values. 





HSE03 


GqWq 


GWo 


GW 


exp. 


£ 


gf.vp. 


SO 


Ge 


0.54 


0.79 


0.82 


0.83 


0.74 


14.0 


16.00 


0.08 


Si 


1.04 


1.32 


1.35 


1.37 


1.17 


9.8 


11.90 




GaAs 


1.12 


1.66 


1.71 


1.75 


1.52 


9.5 


11.10 


0.10 


SiC 


2.03 


2.60 


2.68 


2.76 


2.40 


5.6 


6.52 




CdS 


1.97 


2.55 


2.65 


2.80 


2.42 


4.6 


5.30 


0.02 


AlP 


2.09 


2.69 


2.77 


2.86 


2.45 


6.3 


7.54 




GaN 


2.65 


3.29 


3.38 


3.53 


3.20 


4.6 


5.30 


0.00 


ZnO 


2.11 


2.86 


3.02 


3.33 


3.44 


3.4 


3.74 


0.01 


ZnS 


3.05 


3.69 


3.79 


3.95 


3.91 


4.5 


5.13 


0.02 


C 


5.08 


5.84 


5.92 


6.03 


5.48 


4.9 


5.70 




BN 


5.54 


6.54 


6.66 


6.85 


6.1-6.4 


3.9 


4.50 




MgO 


6.22 


7.94 


8.20 


8.66 


7.83 


2.6 


3.00 




LiF 


11.2 


14.1 


14.5 


15.2 


14.20 


1.8 


1.90 




Ar 


10.1 


13.7 


14.1 


14.7 


14.20 


1.6 


- 




Ne 


14.1 


20.2 


20.7 


21.4 


21.70 


1.2 


- 




MARE 


21 % 


6.8% 


8.0% 


10.0 % 








MRE 


-21 % 


2.3% 


5.3% 


9.3% 











ment of exchange and correlation in the gKS starting func- 
tional improves the prediction of the semi-core d bands, al- 
though a general tendency towards too shallow theoretical d 
states clearly remains. Note, in particular, that the d bands 
shift upwards — away from the experimental values — starting 
from the sX functional. Similar observations have been made 
for an LDAh-U starting point in Refs. 25 and 16. This clearly 
points to a deficiency of the GW approximation. 



B. Selfconslstent QP calculations starting from the HSE03 
functional 



One important question concerns the influence of selfcon- 
sistency and the resulting QP corrections in the case of a gKS 
starting point, especially in comparison to local or semilo- 
cal DFT. For that reason, we have calculated the QP gaps for 
15 materials, without and with partial (only in G) or full (in 
both G and W) selfconsistency with respect to the eigenval- 
ues, starting from the eigenvalues and wave functions of the 
HSE03 functional. In Section IIIIA2I and IIIIA3I this start- 
ing point was found to give the best results for the fundamen- 
tal gaps and the smallest quasiparticle shifts. The materials 
adressed are non-metals, spanning the range from small-gap 
semiconductors to insulators. They are chosen as a subset of 
those considered in Ref. [I3. With the technical details kept 
largely identical, except for the different starting point, the 
data collected in Table |IV] allows for a direct and unbiased 
comparison of the GGA(PBE) starting point"" and the HSE03 
starting point used in this work. The only important differ- 



ence to Ref. 4a is that we now restore the all-electron charge 
density exactly on the plane wave grid for the calculation of 
the correlation energy. This yields technically more accurate 
£/-band binding energies. Details of the applied procedure will 
be published elsewhere. ^^ 

Our results show that the HSE03 gaps, even though they are 
generally closer to the experiment than the DFT-LDA/GGA 
ones, still underestimate the experimental gaps on average 
by 21 %. This underestimation is cured and turned into 
a slight overestimation of about 2.3 % upon the inclusion 
of GqWo quasiparticle corrections. The mean absolute rela- 
tive error (MARE) is reduced to 6.8 %, which is a signif- 
icant improvement compared to the 9.9 % MARE obtained 
for the GGAh-GoWo gaps.'*' Basically, this improvement re- 
sults from the good performance of the HSE03 starting point 
for materials that comprise d electrons such as GaAs, CdS, 
GaN, ZnO, and ZnS, for which the HSE03h-GoWo gap-MARE 
is calculated to be 7.9 %, while it is about 19.2 % in the 
GGA+GqWq approach. We attribute the better agreement to 
the improved description of the pd repulsion on the HSE03 
level, which impacts the energy levels and wave functions (cf. 
Sec. nil A lb . For the rest of the materials, both approaches 
GGA/HSE03 perform on par with a MARE of 4.8/6.1 %. 
While HSE03h-GoWo usually slightly overestimates the band 
gaps, the band gaps for ZnO, ZnS, LiF, Ar, and Ne remain 
underestimated compared to experiment. This may be related 
to the large errors of the HSE03 gap exceeding 30 % for these 
materials. Thus, the inaccuracy of perturbation theory prevails 
for these systems, which share a relatively weakly screened 
exchange (static dielectric constant smaller than 4). Another 
interesting point to note is the different performance of both 
approaches for indirect semiconductors such as Si, SiC, AlP, 
C, and BN. For them the GGAh-GqWo approach performs un- 
expectedly well resulting in a MARE of only 2.6 % while the 
performance of the HSE03 starting point is a little worse on 
average with a MARE of 8.4 %. 

Partial selfconsistency following the GWq scheme is found 
to increase the gaps further by about 0.1-0.2 eV (0.4 eV 
for LiF, Ar, and Ne). Since the HSE03h-GoWo gaps already 
showed the tendency to overestimate the experimental val- 
ues, the MARE increases to 8.0 %. This is in contrast to the 
findings for the GGAh-GWq approach, where an update of the 
eigenvalues in G reduces the MARE to 5.7 %. '^ Following the 
arguments given in Ref. 16, this can be related to the electronic 
macroscopic dielectric constants £„ calculated within RPA, 
which are indicative for the screening involved in Wo. Obvi- 
ously the goo calculated using the HSE03 functional (cf. Table 
lIVI l are underestimated with respect to the experiment and, 
furthermore, they are generally lower than those calculated 
using the GGA XC-functional, due to the in comparison to the 
GGA increased gaps. Only inclusion of excitonic effects, i.e. 
electron-hole binding, allows for the calculation of accurate 
electronic dielectric constants for the HSE03 functionali^^iZ^ 

Further selfconsistent calculations, updating the eigenval- 
ues in both G and W, according to the GW scheme, were per- 
formed. Parallel to the findings for a GGA starting point the 
GW gaps are larger than the GWq gaps, also for the HSE03 
starting point. Consequently they overestimate the experimen- 



TABLE V: Results for the rf-band binding energies of GaAs, GaN, 
ZnO, and ZnS on different levels of quasiparticle selfconsistency 
compared to experimental values starting from HSE03 and PBE 
wavefunctions and eigenvalues, respectively. 





HSE03 


GoWq 


GWo 


GW 


Exp. 


GaAs 


17.2 


17.5 


17.6 


17.6 


18.9 


GaN 


15.4 


16.1 


16.3 


16.5 


17.0 


ZnO 


5.7 


6.1 


6.3 


64 


7.5-8.8 


ZnS 


7.5 


7.2 


7.2 


7.3 


9.0 




PBE 


GoWq 


GWo 




Exp. 


GaAs 


14.8 


16.8 


17.2 




18.9 


GaN 


13.3 


15.4 


16.1 




17.0 


ZnO 


5.2 


6.1 


6.4 




7.5-8.8 


ZnS 


6.1 


6.8 


7.2 




9.0 



tal gaps with a MARE of 10 %, which exceed the GGA+GW 
MARE of 6.1 %. The continued gap increase is found to be 
due to a further reduction of the screening upon updating the 
eigenvalues in W. The difference between the GGAh-GW and 
HSE03H-GW schemes (updating eigenvalues in both G and W) 
can be only related to different starting wave functions. Ob- 
viously the wave functions influence the resulting gap and in- 
crease it on average by 4 %. This is a fairly small change 
confirming the common conjecture that wave functions have 
only a small effect on the band gaps. Similar changes were 
observed in the self-consistent quasiparticle GW (scQPGW) 
suggested by Faleev and Schilfgaarde.i22i22i^ 

Finally, we address the li-band binding energies calculated 
upon the HSE03 and PBE starting point for different levels of 
quasiparticle selfconsistency, as shown in Table IVl We have 
to note that the present calculations do not include the s and p 
orbitals with the same main quantum number as the semi-core 
d shell. However, since here the core-valence interaction is 
approximated by HE exchange rather than by LDA, as inher- 
ent to conventional pseudopotential calculations, the values 
given in TablelVJalreadv provide a reasonable estimate for the 
(i-band binding energies (cf. Ref. 18, and 16). Updated val- 
ues for the PBE case are also supplied; the present values su- 
percede those in Ref. 16 and are more accurate, since the all- 
electron charge density is now accurately restored on the plane 
wave grid.^^ In general, the QP binding energies increase over 
the HSE03 one-electron values due to the quasiparticle cor- 
rections, which can be understood mainly from the effects of 
the enhanced core-valence interaction in the GW calculations. 
Only in the case of ZnS this trend does not hold for yet un- 
known reasons. However, compared to the experimental val- 
ues all three quasiparticle schemes studied here underestimate 
the (i-band binding energies. With increasing selfconsistency 
along the row (GqWo, GWq, and GW) the d bands shift to 
larger binding energies. Thereby, for the Ga 3d levels, the 
calculated values approach the experimental ones. In the case 
of the Zn 3d levels, which are just below the p-like upper va- 
lence band complex, the situation is different with a more pro- 
nounced underestimation of the binding energies. Analogue 
observations have been made for a GGA starting point, ^^'^^i^ 
which gives a smaller d-band binding energy than the HSE03 



starting point for each level of selfconsistency. We note that 
Fleszar and Hanke observed that the inclusion of vertex cor- 
rections in the self-energy shifts the d states to stronger bind- 
ing energies, "" suggesting that the neglect of such corrections 
is responsible for the erroneous behavior of the GW approxi- 
mation for d states. 



IV. SUMMARY AND CONCLUSIONS 

We have presented GqWq QP calculations starting from a 
variety of XC functional: LDA, sX, HSE03, PBEO, and HF 
for Si, InN and ZnO. We have shown that the gKS schemes, 
which take into account a screened exchange potential or part 
of it (sX, HSE03), give rise to eigenvalues close to the QP ex- 
citation energies. The resulting GqWo corrections were found 
to yield QP energies in good agreement with the experimental 
data, and the QP shifts are less dispersive than those calcu- 
lated upon LDA. Overall the HSE03 and PBEO functional 
gave one-electron energies very close to the successive GW 
calculations, resulting in small QP gap corrections across the 
considered energy range. For the HSE03h-GoWo case, the final 
QP energies were in very good agreement with experiment, 
whereas the PBEO functional was found to yield too large QP 
gaps. We traced this back to a significant underestimation of 
the screening for the PBEO functional, when the random phase 
approximation is used (also applied to determine W). 

Furthermore, the QP gaps for 15 materials comprising 
small and large gap systems were calculated, in order to pro- 
vide a benchmark of the HSE03 starting point against the 
GGA one. It was shown that the HSE03h-GoWo approach 
yields an almost halved overall error for the fundamental gaps 
compared to the GGA starting point. The largest improvement 
over the LDA/GGA starting points were found for materials 
with shallow d states such as ZnO, ZnS, InN, GaAs, and GaN, 
where the LDA/GGA starting point suffers from a significant 
underestimation of the tZ-band binding energies and a conse- 
quently overestimated repulsion between p- and d-like states. 
Since the tZ-band binding energies calculated using one of the 
gKS schemes (e.g. HSE03) are closer to the experiment, the 
influence of the pd repulsion on the gap is described more 
accurately. 

Furthermore, the effects of different degrees of selfconsis- 
tency were investigated. It was found that both, selfconsis- 
tency in G (GWq) and selfconsistency in G and W (GW), im- 
pair the agreement with experimental data. For GWq, this is 
in contrast to the findings for a GGA starting point. ^^ This 
could be traced back to the poorer description of the dielec- 
tric screening for the HSE03 starting point. Selfconsistency 
according to the GW scheme with an update of eigenvalues 
in both G and W further diminishes the agreement with mea- 
surements, due to a further reduction of the already underes- 
timated screening. This is analoge to the findings for GGA 
based GW calculations. In general, concerning selconsistent 
QP schemes, the present work confirms the observation al- 
ready made in previous work: to obtain accurate QP gaps it is 
essential to use an electronic response function and a screened 
interaction W that agree closely with experiment, ^^'^^ and to 



combine this screened interaction with an accurate Green's 
function G. For the GWq case based upon HSE03 wave func- 
tions, the overestimation of the gaps clearly relates to an 
underestimation of the static electronic screening employing 
HSE03 and the random phase approximation. In this light, 
the success of the perturbative single-shot HSE03h-GoW(3 ap- 
proach is a little bit fortuitous, since the overestimation of 
the screening is partially canceled by too small gap correc- 
tions obtained using the single-shot perturbative GqWq ap- 
proach. Probably the same is true for some other single shot 
approaches, such as EXX-OEPh-GoWq. Nevertheless, if com- 
putational efficiency is an important issue — and it more often 
is than not — then the HSE03h-GoW() approach is indeed an 
excellent balance between accuracy and speed. The calcula- 
tions are as efficient as for the commonly used LDAh-GoWq 
method, and, with very few exceptions, the errors are smaller 
than 10 %. If better accuracy is required, it can be achieved, 
but only by updating the wave functions and including exci- 
tonic effects in the calculation of the screening properties, i.e., 
vertex corrections in WP For most mater, such calculations 
are currently not feasible due to the large computational re- 
quirements. 

A similar accuracy as for the HSEh-GqWo approach can be 
achieved by starting from GGA wave functions and eigen- 
values, and updating the eigenvalues in G until convergence 
is reached (see Ref. 16). This approach yields comparable 
errors as HSE03h-GoWo, but it is computationally more de- 
manding, since several iterations are required to converge G. 
The latter method is also problematic for materials with an in- 
verted band order in LDA/GGA. Which approach to choose 
(HSE03h-GoWo or GGAh-GWo) is to some extend a matter of 
taste, and the final results are usually very close and often 
bracket the experiment. If efficiency and robustness (band 
order) are issues, the HSE03h-G()Wo approach seems to be 
preferable, and it is certainly much more accurate than the 
traditional LDAh-GqWo method. 

Concerning the position of the d levels, the HSE03h-GoWo 
method shows an underestimation of the cZ-band binding ener- 
gies by about 1 eV for almost all materials. Similar observa- 
tions were made for the LDA/GGA case ^'V^i^' and LDAh-U 
based GW calculations, '^-^^ and self-consistent quasiparticle 
GW (scQPGW) calculations.^^ The underestimation of the d- 
band binding energy is thus universal to the GW approxima- 
tion and not related to the starting wave functions. The origin 
for this underestimation is yet unknown, but it is most likely 
related to the fact that the GW approximation is not entirely 
free of self-interaction errors, and only inclusion of vertex cor- 
rections in the self-energy might remedy this deficiencyi^ 
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